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We investigate the question of stability of a solid thin film which experiences external inter- 
actions such as van der Waals forces from a contacting surface or forces from an external 
electric field. Both perfectly elastic and viscoelastic material behaviours are considered in 
linear stability analysis performed here. These analyses indicate that for sufficiently soft 
(shear modulus between 1 and 10 MPa) and nearly incompressible films (Poisson's ratio 
close to 0.5), bifurcations are possible, i. e., the surface of the film becomes non-planar. 
The modes of bifurcation and rates of growth of perturbations are determined as a function 
of material parameters. The results of this study are of significance in understanding the 
adhesive properties between a soft material (such as rubber) and a comparatively rigid solid 
(such as steel), and the behaviour of soft solid films in an electric field. 



Instabilities and pattern formation in thin solid and 
liquid films are of interest both from a scientific and 
technological view point. Morphological instabilities 
in thin liquid films occur due to causes such as com- 
petition between capillary forces and van der Waals 
interactions |jj or an external electric field Q and can 
often lead to dewetting leading to interesting patterns. 
Morphological instabilities are also common in solid 
films; for example, in stressed solid films the strain en- 
ergy drives the surface roughening in competition with 
the surface energy with surface mass diffusion being 
the dissipative mechanism || ||. 

Analysis of interacting thin films has hitherto been 
restricted to fluid films. Here, we pose the question of 
stability of a thin solid film bonded to a rigid substrate 
whose free surface experiences an effective force. This 
force may arise from any of the various causes such as 
a van der Waals interaction with another contacting 
surface nearby and/or with the substrate, an external 
electric field, etc. The theoretical analysis presented 
in this paper indicates that for a soft and nearly in- 
compressible solid thin film, instabilities are possible 
and that the film "buckles". Physically, this instabil- 
ity occurs because it is possible, for sufficiently large 
interaction forces, to reduce the net potential energy 
of the system (the elastic strain energy and the surface 
energy of the film + potential energy of interaction of 
the surface) by a periodic non homogeneous deforma- 
tion in the film. We believe that these results could 
be useful in understanding phenomena of adhesion be- 
tween materials (such as rubber and steel), behaviour 
of thin films in an electric field etc. 

The system considered here is shown in fig. [j] - a 
film of height h bonded to a rigid substrate described 
by coordinates {x\,X2) such that surface of the film 
S interacting with external agency has x-i — and 
that bonded with the rigid substrate has X2 = —h. 
We restrict attention to plane strain deformations of 




FIG. 1: Film bonded to a rigid substrate. The surface 
of the film experiences external forces. 

the film for the sake of mathematical simplicity and to 
understand the essential physics. The total potential 
energy of this system system is 

J W{e)dV + J ( 7v /l + (M 2 ,i) 2 -(7(u-n)) dS (1) 

where e is the strain tensor, W(e) is the elastic strain 
energy density, 7 is the surface energy, U(u ■ n) is the 
interaction potential between the surface of the film 
and the external agency such as a contactor or an elec- 
tric field, u is the displacement vector and n is the 
outward normal to the surface. Linearised analysis is 
performed by expanding the interaction term U(u ■ n) 
in a power series about u = and retaining all terms 
up to quadratic order in u. The resulting approximate 
energy functional is 

/ W(e)dV + / 71/I + («2,i) 2 dS 
Jv J s 

- f (u + F o u-n + ^Y(u-nf^dS (2) 
J s ~ 

where 

U = U(0), F = U'(0) and Y = U"(0). (3) 
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FIG. 2: Bifurcations modes (hk) as a function of hY/ fi 
for various values of v with j/(J>h = 0. (b)Bifurcations 
modes (hk) as a function of hY/fj, for various values of 
/ fih with ^ = 0.4. 

The equilibrium stress field cr in the film (which min- 
imises the potential energy (^) over an appropriate 
length of the film) satisfies the equilibrium equation 
V • cr = in V and the boundary condition 



a ■ n = jU2,un + F n + Y(u ■ n)n 



(4) 



on S. Taking the film to be an isotropic linear elastic 
solid with shear modulus \i and Poisson's ratio v, gives 
a standard expression for the strain energy density |B| 
with a resulting expression for the stress tensor ex- 
pressed in terms of the gradient of displacement. Thus 
the problem can be cast into a boundary value problem 
for the unknown displacement field with the boundary 
condition of vanishing displacements at 22 = — h at the 
the film substrate interface in addition to (Q). 
The Homogeneous Solution: A solution to the above 
boundary value problem exists such that the stresses 
in the film are equal everywhere. This homogeneous 
solution (u h ) is u\ = everywhere, and u\ has a linear 
variation with X2 starting from at X2 — —h, i.e., 



u 2 (xi,Q) — u 



For the case when v 



Fa 



- Y 



(5) 



( 2(l-i/) M 
I (l-2v)h 

0.5, i. e., the incompressible 



limit, the homogeneous solution is such that the dis- 
placement vanishes everywhere in the film, and a pres- 
sure field p develops such that p(x±, £2) = F . So long 

as 



Y < Y m 



hY m 



2(1-^) 
(1 - 2v) 



(fi) 



the homogeneous solution is meaningful in that u has 
the same sign as F . This conditions on Y is most easily 
met when v is close to 0.5 (the r.h.s. of (||) tends to 00 
as v tends to 0.5), i. e., when the material in nearly 
incompressible. It is this class of materials that the 
focus of this paper. Nevertheless, results are presented 
for all values of v for the sake of completeness. 
Bifurcations: What are the conditions (on Y, fj,, v, h, 7) 
for another solution (inhomogeneous state) to exist? If 
such a solution exists, it can be taken to be of the form 
u h + u, where the symbol u now stands for a "bifur- 
cation" displacement field. This bifurcation field must 
satisfy the equilibrium equations in the bulk and the 
rigid boundary condition at the film substrate inter- 
face, just as the homogeneous solution. On the surface 
of the film at X2 = 0, the bifurcation field satisfies (here 
cr is the additional stresses due to u), 



an — 7112, iin + Y(u ■ n)n, 



(") 



instead of (Q). To investigate the existence of a non- 
trivial solution to the problem defined above, the bi- 
furcation fields are assumed to have the form 



Uj(xi,X2) = e lkxi Uj(x2) 



(8) 



where k is a real positive wavenumber. The problem 
of finding nontrivial bifurcation fields can be cast into 
the problem of finding those values of k such that the 
functions Uj(x2) are nontrivial. It can be shown that (a 
detailed account will be published elsewhere) nontrivial 
bifurcation fields of the form (Ji) exist for those values 
of k that satisfy the equation 

(k [4e 2hh hk 2 (V - (1 - v)i) + (e ihk - l)fc 7 (3 - 7u + Au 2 ) 

+a ((3 - 4 M )(1 + e 4hk ) - 2e 2hk (5 - I2v + 8u 2 ))] ) j 

((1 - v) [(3 - 4v)(e 4hk - 1) - 4hke 2hk ] \ = Y (9) 

This relation is valid for the incompressible case as 
well (i. e., when v = 0.5). Real roots of (M) are sought 
when Y < Y m which is the range of Y for which the 
homogeneous solution is valid. 

We first focus attention on the case when 7 vanishes. 
Fig. ||a depicts graphically the solution to @, i. e., for 
a given value of v, the values of k that solve (||) are 
plotted as a function of Y (hY/ fi in non-dimensional 
terms) . The important results may be noted: (i) There 
are no bifurcation modes for any value of v when 
hY/fj, < 2. (ii) For all values of v, k — is a bi- 
furcation mode when Y = Y m . (iii)When v < 0.25, 
there are no bifurcation modes for Y < Y m . (iv) When 
v > 0.25, there are two modes starting from a critical 



3 



value Y c (such as the point C shown in fig. ||a) that 
depends on the value of v until Y reaches Y m . When 
the film is incompressible hY c /fj, = 6.22 and the cor- 
responding bifurcation mode has hk c = 2.12. For this 
case bifurcations are possible for all values of Y greater 
than 6.22fi/h, with two possible values of k as shown 
in the fig. |^a. 

Next, we consider the case when 7 / 0. Fig. |^b 
shows a plot of the possible wavenumbers of bifurcation 
modes for various values of 7 with v = 0.4. The key 
effect of the surface energy on the bifurcation modes 
are noted as follows: (i) Surface energy inhibits bifur- 
cation, in that a larger value of Y c is effected with a non 
zero value of 7. The critical mode k c decreases with in- 
creasing 7. Both of these results are as expected since 
a larger value of k implies a larger energy penalty in 
terms of surface energy, (ii) As 7 gets larger Y c ap- 
proaches Y m . In fact, it can be shown that Y c equals 
Y m when 7 = j m where 

2u(iu - 1) 



2.2 



7m 



3(1 -2v) 



(10) 



a result which is pertinent when v > 0.25. The curve 
for 7//1/1 = 4.0 for the case of v — 0.4 shown in fig. 0, 
graphically illustrates this point. If 7 > j m , then there 
are no bifurcations in the physically meaningful range 
Y <Y m . 

A more detailed analysis gives the following formulae 
for Y c and k c as a function of 7 and v when 7/ iih <C 1 
and v — > 0.5: 



-Yc(v,f/nh) 



■ .22 - 10.46(1 - 2t 
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hk c (u,j/fj,h) = 2.12 - 2.86(1 - 2v) - 2A2-. (11) 
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FIG. 3: Bifurcations modes (hk) as a function of hY/fi 
for various values of v with 7////1 = 0. (b) Bifurcations 
modes (hk) as a function of hY/ fx for various values of 
j/fth with v = 0.4. 



It is also interesting to consider the time evolution of 
deformation in the film so as to obtain the dominant or 
the fastest growing mode. To this end, the film is con- 
sidered to be viscoelastic with a constitutive relation 
of the form 



Vu T ) + 



-V ul 



2 M Q(V U 
+2n Q(Vu + Vu T ) - iv ■ ul 



(12) 



where (') stands for the time derivative, 77 is a viscosity 
parameter and I is the second order identity tensor. 
In the consideration of the time evolution of the sys- 
tem, inertial effects are neglected since the time scale 
of interest is much larger than the time scale of the 
propagation of an elastic wave through the thickness 
of the film. 

The Homogeneous Viscoelastic Solution: The homoge- 
neous solution of the field equations with the viscoelas- 
tic constitutive relation (IT4) is 

ul = 0, u%(xi,X2,t) = u (l + j^j (l — e w *^ (13) 

where u h is given by 

h __3_ / 2(1 - v)n 
U An \(l-2v)h 



Y 



-^(Y m -Y). (14) 



From (14) it is evident that the time dependent ho- 
mogeneous solution tends to the elastic homogeneous 
solution as t — ► 00 when Y < Y m . If Y > Y m , the 
present analysis indicates that the homogeneous solu- 
tion blows up as t — > 00. 

Growth of Perturbations: Just as in the case of the 
elastic film, it is of interest to investigate the growth of 
perturbations of the homogeneous solution. The per- 
turbations u are assumed to be of the form 



Uj(x\,X2,t) = e xl Uj(x2)e u 



(15) 



For a given k, the rate of growth ui is determined by 
insisting that the the perturbation satisfies equilibrium 
equations and boundary conditions and that they be 
nontrivial. The relation between u and k can be ob- 
tained by replacing /1 and v in (|^) respectively by fx* 
and v* where, 



u = fj, + nui, 



— (1 — 2v)r]0J 
3/it + (1 — 2u)nuj 



(16) 



This procedure results in a cubic equation for to. 
The solution of this equation is obtained by numeri- 
cal means. 

The solution for uj indicates that for Y c < Y < Y m) 
all perturbation modes with wavenumbers between the 
two bifurcation modes given by the elastic analysis 
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FIG. 4: The stable and unstable regions in the parame- 
ter space. Region /: Homogeneous solution stable. Re- 
gion II: Homogeneous solution unstable with two elas- 
tic bifurcation modes and a fastest growing viscoelastic 
perturbation. Region 777: Homogeneous solution "un- 
physical" . 

are unstable i. e., u for these modes are positive. In- 
deed, there is a mode with wavenumber (k m ) between 
wavenumbers of the two elastic bifurcation modes such 
that the rate of growth (uj) is a maximum. Fig. |^(a) 
shows a plot of k m as a function of Y (Y c < Y < Y m ) 
for various values of v (with j/f^h — 0). When v < 0.5, 
the value of k m starts at k c when Y = Y c and monoton- 
ically falls with increasing Y. For the case of v = 0.5, 
k m = k c for all values of Y. When 7 ^ 0, k m is smaller 
as is evident from fig. ^(b) ; the effect of surface energy 
on the fastest growing mode becomes increasingly less 
significant for large values of Y. Just as in (111), an 
analytic result can be derived for k m for small values 
of y/fih, v -> 0.5 and h(Y - Y c )/fx < 1: 



hk m (u, — ) = hk c (u. 
fin 



Hh J 



+ [ 0.39^: - 0.46(1 - 2v) ) -(Y - Y c ) (17) 



Instability in a thin film whose surface experiences 
forces depends on three key sets of non-dimensional pa- 
rameters namely the Poisson's ratio v : the normalised 
second derivative of the interaction potential hY/ [i and 
the normalised surface energy "f/fJ-h. The whole pic- 
ture of stability and bifurcation in this system and its 
dependence on the nondimensional parameters can be 
depicted pictorially as shown in fig. [I| Region I in fig. [| 
is where the homogeneous solution is unique and stable 
while region marked III in the figure corresponds to 
the case when the homogeneous solution is "unphysi- 
cal", i.e., this analysis is not adequate. Region II is 
the most interesting - this corresponds to nearly in- 
compressible material behaviour. In this region the 
homogeneous solution is unstable, with two possible 
elastic bifurcation modes; a viscoelastic analysis pre- 
dicts a fastest growing mode with a wave vector that 
lies between the two elastic bifurcation modes. 

We now turn to specific cases of the type of sys- 



tem considered in this paper. First, we consider a 
rigid contactor interacting with the film via van der 
Waals forces. Assuming that the contactor is at a 
distance d above the undeformed surface of the film, 
the interaction potential U can can be taken to be 

U ( u ■ n ) = i2*(u A n-d)> with F ° = 53*> Y = 2^F- 
Taking the film to be made of rubber (fi = 1 MPa, 
v = 0.5, 7 = 0.1J/m 2 ) and h = lmicron with A w leV. 
When d = lOnanometers we get hY/ fi = 1.6 and for 
d = 5nanometers, hY/ /i = 25.6. Since the latter value 
is greater than hY c /fi which is 6.63 when j/ixh = 0.1 
(which is the present case), it is clear that the con- 
dition for bifurcation will be achieved as d is reduced 
from lOnm to 5nm. Thus as the contactor approaches 
the film, the film would buckle. This implies that the 
contact that forms between the contactor surface and 
the film will not be planar. We are not aware of any 
experimental work that can corroborate our results. 
We do, however, hope that the contents of this paper 
will be useful in designing experiments to verify our 
conclusions. 

The second case considered is that of a film interact- 
ing with an external electric field. The system consists 
of two plates separated by a distance d] the bottom 
plate is coated with a nearly incompressible polymeric 
film of height h. A potential difference of V is applied 
between the two plates. The quantity of interest is the 
value of the gap thickness d — h at which instability 
occurs in the film. The potential of interaction for this 

case is given by U{u ■ n) = 2(epd _ ( E e p l P 1 ^ h+u . n)) where 
£0 is the permittivity of free space, e p is the dielectric 
constant of the polymer. Taking the mechanical prop- 
erties of the polymer to be same as in the previous case, 
and taking e p — 3, we get that the critical gap thickness 
d — hoi 0.05micron for a film of height 0. lmicron with 
the applied voltage of 100V. A gap thickness smaller 
than 0.05micron will cause the film to buckle. It is ev- 
ident that large electric fields are required to cause the 
instability. 

We do wish to point out that this analysis is based 
on a linearised model, and will only provide the modes 
of instability, i. e., the wavelength of surface undula- 
tion and not the magnitude. A nonlinear analysis is 
required to obtain such a quantity and will be pursued 
in subsequent papers. 
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